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Distributions of eigenmodes are widely concerned in both bounded and open systems. In the 
realm of chaos, counting resonances can characterize the underlying dynamics (regular vs. chaotic), 
and is often instrumental to identify classical-to-quantum correspondence. Here, we study, both 
theoretically and experimentally, the statistics of chaotic resonances in an optical microcavity with 
a mixed phase space of both regular and chaotic dynamics. Information on the number of chaotic 
modes is extracted by counting regular modes, which couple to the former via dynamical tunneling. 

The experimental data are in agreement with a known semiclassical prediction for the dependence of 
the number of chaotic resonances on the number of open channels, while they deviate significantly 
from a purely random-matrix-theory-based treatment, in general. We ascribe this result to the 
ballistic decay of the rays, which occurs within Ehrenfest time, and importantly, within the timescale 
of transient chaos. The present approach may provide a general tool for the statistical analysis of 
chaotic resonances in open systems. 


Introduction. The statistics of chaotic resonances has 
been a central topic of theoretical and experimental in¬ 
terest for decades [1-8] , as a doorway to better under¬ 
stand chaotic scattering in quantum mechanics [9-12]. 
Counting chaotic resonances notably finds applications 
in optical resonators [13-15], where chaos can be used, 
for example, to enhance energy storage [16] or enhance 
coupling [17, 18]. In the realm of fundamental prob¬ 
lems, an estimate of the number of states within a certain 
frequency interval is given by the fractal Weyl law [19- 
23]. Predicted more than a decade ago, this theory is 
still awaiting experimental confirmation at optical fre¬ 
quencies. The delay is mainly due to two obstacles: (i) 
the theoretical problem of accounting for partial absorp¬ 
tion [24, 25]; and (ii) the experimental challenge of ana¬ 
lyzing overlapping resonances [26] . By introducing meth¬ 
ods to overcome these hurdles, we theoretically and ex¬ 
perimentally study the statistics of chaotic resonances in 
two-dimensional asymmetric optical microcavities, and 
thus present a significant result in the context of fractal 
Weyl laws at optical frequencies. An absorber is placed 
underneath the microcavity appositely to realize a vir¬ 
tually full opening, and therefore solve the problem of 
partial absorption. In order to handle overlapping res¬ 
onances, we exploit the mixed classical phase space of 
the present cavity, and exclusively count high-Q regular 
modes, easily recognized in the measured spectra due to 
their narrow linewidths, and that enables us to draw in¬ 
formation on the low-Q chaotic modes, coupled to the 
regular ones via dynamical tunneling [27]. 

Theoretical model. The dynamics inside the cav¬ 
ity is described by a non-Hermitian Hamiltonian H = 
Hq + V, where Hq has eigenstates \C^) (regular) and 
\Cn) (chaotic), while V represents the coupling be- 
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FIG. 1. (color online) (a) Sketch of the deformed microcav¬ 
ity with an inner absorber, characterized by the angle Oa- 
(b) Poincare surface of section of the microcavity (0 = 
with deformation factor 77 = 11.7%. Blue orbits at the top are 
regular, black ones are chaotic. The red solid line indicates 
the angle of total internal reflection, while the dashed curve 
is given by an absorption angle such that r ~ 0.77. Different 
shades of color indicate loss to the absorber (lighter) and by 
refraction into air (darker), (c) Schematic representation of 
the free-space coupled cavity system, (d) A typical transmis¬ 
sion spectrum with the high-Q regular modes highlighted. 


tween them [28, 29]. Following a standard approach 
[30, 31], the electromagnetic field excited by the incident 
beam can be written as a superposition of one regular 
and several chaotic modes, fj{x,t) = a^(t)c^(x)e“*‘^^ + 
^^bn{t)cn{x)e~^^^^. The a^j and bn are oscillator am¬ 
plitudes driven by the laser beam, coupled to each other 
with strength Vn (real), with damping rates juj and 7n, 
respectively. Then, in the overdamped regime (|Yn| ^ 
7n), the dynamical equations for the time-dependent en- 








velopes take the form [17] 
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where fn is the coupling strength of the n-th chaotic 
mode with the laser beam of amplitude Eq and frequency 
uq. We are interested in the steady-state solution, ob¬ 
tained by setting = 6 ^ = 0. The amplitude a^j of the 
envelope of the regular mode is then found to be 

EoEfn^ 

aoj = - - -( 2 ) 
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We only sum over chaotic modes with linewidth 
smaller than a certain value set by 7 , and use the averages 
/, y ,7 to approximate the summations as 

[32]. The result is, at resonance 
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where e = Eof/V, and T = • 

Equation (16) is central to the present work, as it 
bridges between the number of excited regular modes 
(proportional to la^^p), which we measure directly, and 
the number of chaotic modes which we estimate. In 
the spirit of Weyl law, we test the prediction [33] 
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for the dependence of on the number of open chan¬ 
nels N . Here M is the total number of chaotic states, 
Td = M/N is the mean dwelling time of a ray in the cav¬ 
ity, and yU is of the order of the Lyapunov exponent of the 
chaotic dynamics [32] . We focus on the role of the prefac¬ 
tor which accounts for the instantaneous-decay 

modes escaping from the system within the Ehrenfest 
time of quantum-to-classical correspondence [34] . By re¬ 
moving the prefactor, the remainder of Eq. (4) is solely 
based on random matrix theory (RMT) [35], and will also 
be tested against the observations. 

Absorber in the optieal eavity. In order to achieve the 
full opening required to test these predictions, we intro¬ 
duce an absorber in the cavity. In the analysis, the di¬ 
electric microcavity [Eig. 1(a)] has the deformed circle 
p{(j)) as boundary [32], and encloses an absorber of shape 
p{(f)) — R. Eigure 1(b) shows the classical phase space, 
together with the critical line of total internal reflection 
(sin 6 >c), as well as the line given by the incidence angle 6 >a, 
below which the reflected ray hits the absorber [36] . We 
assume that the rays hitting the absorber are completely 
absorbed. The rays that escape the cavity by refraction 
into the air with an angle of incidence 0 Oc are very 
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FIG. 2. (color online) (a), (b) Number of chaotic states vs. 
the rescaled absorber-to-cavity ratio obtained by RMT and 
Weyl law, respectively, (c), (d) The corresponding expecta¬ 
tions for |acu|^. In (d), the smaller /x, the more visible the tail 
of the excitation probability. The red vertical line corresponds 
to the critical angle, sin^c — 0.69. Here sin^th = 0.6. 


lossy, and, as such, they are not expected to contribute 
to the excitation of the regular modes, consistently with 
Eqs. (14) and (16). 

Eor that reason, we only take into account the states 
supported on a strip of the chaotic phase space with mo¬ 
mentum above a certain threshold, smO > sin^th, to 
be chosen below but close enough to the critical line of 
total internal reflection. Let ns introduce the notation 
^ = sin Oa — sin ^th to indicate the strip of the phase space 
opened by the absorber. Using the picture of Ref. [33], 
there are effectively N open channels out of the M Planck 
cells available in the phase space, produced by the ab¬ 
sorber (full opening, Na) and the refraction out of the 
cavity (partial opening, W), so that the mean dwelling 
time of a ray is given by 
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Td = 


NaENr' 


( 5 ) 


with Nr = ^ d sin 0T{sm0), T transmission coef¬ 

ficient according to Eresnel law, and A area of the phase 
space in exam, while Na = M^jA. 

Recalling the original purpose of studying the statis¬ 
tics of chaotic resonances, we proceed by steps and 
first examine the RMT-based prediction Ur^^RUT — 
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, rewritten as (set ^ = l/rd = + 
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The theoretical expectation is portrayed in Eig. 2(a): 
^ 7 ,RMT changes relatively slowly for ^ small such that 
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when the loss is mainly due to refraction into 
air. Otherwise decreases more rapidly and lin¬ 

early with ^ in the region of total internal reflection, when 
the loss is entirely due to the absorber. Plugging Eq. ( 6 ) 
into Eq. (16), the probability of excitation of the high-Q 
regular modes I^^wIrmt starts to fall off as ^ reaches some 
critical value, controlled by the parameter 7 [Eig. 2(c)]. 
The other parameter E = 7 ^ 7 /ME^ controls the slope of 
the curve. It is noted that the probability of excitation of 
the regular modes decreases dramatically when the loss 
is entirely due to the absorber, in which case the system 
is fully open. 

On the other hand, the semiclassical estimate (4) be¬ 
comes, as a function of 


%,Weyl 





i 

1 - e-'f 


( 7 ) 


The quantity /i, of the order of the Lyapunov exponent 
of the chaotic region of the phase space [32], is what 
really characterizes (7), which resembles the linear RMT 
prediction ( 6 ) for large enough /r, and otherwise becomes 
visibly nonlinear [Eigs. 2(a) and 2(b)] when /i ^ 1 . This 
nonlinearity produces a characteristic tail in the curve 
expressed by Eq. (16) [Eig. 2(d)], meaning that the effect 
of the Ehrenfest time scale on the excitation of the regular 
modes is most evident slightly above the onset of chaos. 

Experimental setup and measurement. The experimen¬ 
tal apparatus consists of a deformed toroidal microcavity 
[18] coupled to a laser beam (of wavelength A 1550 or 
635 nm), as shown in Eig. 1(c). The microtoroid (refrac¬ 
tive index 1.44,1.46 depending on A) has principal (mi¬ 
nor) diameters of 120 /im(5 /im), consistently with the 
two-dimensional model [32]. Thus the effective Planck 
constant heff ^ A/a ^ 10“^ (a: principal diameter) jus¬ 
tifies the semiclassical analysis. The microcavity is fabri¬ 
cated through optical lithography, buffered HE wet etch¬ 
ing, XeE 2 gas etching, and CO 2 pulse laser irradiation. 
The resulting silica microtoroid is supported by a sili¬ 
con pillar of similar shape, which has a high refractive 
index 3.48,3.88), and it acts as the absorber in the 
model. After each measurement of the free-space trans¬ 
mission spectrum [Eig. 1(d)], the top diameter of the sili¬ 
con pillar, connected with the silica disk, is progressively 
reduced by a new isotropic XeE 2 dry etching process. 
In this way we control the openness of the microcavity 
with the ratio r between the top diameters of pillar and 
toroid. Einite element method simulations show that the 
light power decreases to less than 5% of the input value, 
when propagating by a distance of 20 /im inside the 2 - 
/im-thick silica waveguide bonding with a silicon wafer 
[32], as is reasonable to expect, given the high refrac¬ 
tive index of the silicon. Thus the silicon pillar acts as 
a full absorber, consistently with the present model. On 
the other hand, high-Q regular modes living inside the 
toroidal part, whose cross section has minor diameter of 
5 /im, do not leak into the silicon pillar and therefore are 
not directly affected by the pillar size [32]. The depen- 
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FIG. 3. (color online) Normalized transmission and top-view 
optical images of the cavity with r ~ 0.81 [(a) and (e)], 0.77 
[(b) and (f)], 0.70 [(c) and (g)], and 0.64 [(d) and (h)]. Inset of 
(a) shows background noise. Insets of (b)-(d) show the high- 
Q modes. Reflection of the silica-to-silicon interface results 
in a brighter color for the silicon pillar in the optical image 
(boundary shown by red dashed curves). Scale bar is 50 /xm. 



dence of the free-space transmission spectra on the pillar 
size is shown in Eig. 3. When the pillar approaches the 
inner edge of the toroid [Eigs. 3(a) and 3(e), r 0.81], no 
high-Q regular modes are observed in the spectrum, since 
most of the probe laser field in the cavity radiates into 
the silicon and cannot tunnel to high-Q regular modes. 
As we gradually reduce the size of the pillar [Eig. 3(f), 
r 0.77], increasingly many high-Q modes appear in the 
spectrum [Eig. 3(b)]. When the absorber-to-cavity ratio 
r is small enough [Eigs. 3(g) and 3(h), r < 0.7], the trans¬ 
mission no longer changes sensibly [Eigs. 3(c) and 3(d)], 
and the number of high-Q modes in the spectrum also 
stabilizes. 

Statistics of chaotic resonances. As anticipated, we use 
the transmission spectra to test the theory, by counting 
the excited high-Q regular modes for different sizes of the 
silicon pillar [37]. The results are illustrated in Eig. 4, for 
two microcavities of distinct deformations. In the RMT- 
based approach [Eqs. (16) and ( 6 )] we have two fitting 
parameters, 7 and P, while the total number of chaotic 
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Rescaled absorber-to-cavity ratio ^ 


FIG. 4. (color online) (a), (b) Number of high-Q regular 
modes {uuj) observed experimentally in the microcavities of 
rj = 11.7%, 6% respectively, vs. rescaled absorber-to-cavity 
ratio Blue dashed and red solid curves are respectively 
RMT and Weyl-law best fits, whose goodness is assessed by 
Here sin^c — 0.69 and sin^th = 0.6. Inset: area where 
the two curves differ most. 


states is estimated theoretically as M ^ area 

of the phase space we consider). Figure 4(a) shows over¬ 
all agreement between the experimental data (dots) and 
this theory (blue dashed curve), which however deviates 
from the tail visible at larger sizes of the absorber. The 
discrepancy becomes more apparent for the cavity with a 
lower deformation factor rj [Fig. 4(b)]. A quantitative test 
of the goodness of the fit yields a reduced — 1-7, 5.6 
for T] 0 ^ 11.7%, 6 % respectively [32, 38], the latter devi¬ 
ating significantly from the optimal value of unity. The 
fitted maximum escape rate 7 0.2,0.3 is the inverse 

(in units of Poincare time) of the minimum escape time 
Tesc of the chaotic rays contributing to the excitation of 
the regular modes, from which Q = 27ruresc r\j 10 ^, on 
average (z/ is the frequency of the laser beam). We find 
this estimate consistent with the typical order of Q inde¬ 
pendently obtained from ray-dynamics simulations [32], 
which suggests the fitted parameter makes physical sense. 

Next, we test the semiclassical correction (7), using the 
finite time Lyapunov exponent /i evaluated by direct it¬ 
eration [32, 39], with the estimated parameter M and the 
fitted parameters 7 and F. It is found that the semiclas¬ 
sical correction (red solid curves) fits the experimental 
data better than the RMT-based estimate, especially at 
the smaller deformation, where the two predictions dif¬ 
fer the most due to the smaller /i [cf. Fig. 2 ]. Here 

— 1-3,1, indicating good agreement. In particular, we 
are now able to account for the tail of the curve, which 
corresponds to the microcavity having the largest open¬ 
ings and thus with the maximum number of instanta¬ 
neous decay states, where the semiclassical correction is 
decisive. More experimental and fitting results at dif¬ 
ferent wavelengths and deformations support the above 
explanation [32]. 

Conclusion and discussion. Let us summarize the work 
done and the results obtained. By counting the high-Q 
regular modes excited via dynamical tunneling as a func¬ 
tion of the number of open channels, we have studied the 
statistics of the chaotic resonances in a dielectric micro¬ 
cavity with a full absorber. 



FIG. 5. (color online) Survival probability in the chaotic re¬ 
gion (logarithmic scale). Points: average survival probability 
P{t) of a ray in the microcavity vs. t (in units of Poincare 
time) at ?7 = 11 . 7 %,^ = 0 . 13 , from 10 ® randomly-started 
trajectories. Line: P{t) = exp(—t/r^), ra — 6. Inset: the 
long-time simulation showing algebraic decay. 


As main result, the experimental data deviate from a 
purely RMT-based prediction, while they exhibit better 
agreement with a semiclassical expression that factors 
out the number of instantaneous decay modes. Impor¬ 
tantly, the latter estimate depends on the Lyapunov ex¬ 
ponent of the chaotic dynamics, and it accounts for a 
characteristic tail in the decay of the number of regular 
modes, which we interpret as a signature of the ballis¬ 
tic escape of the rays into the absorber, occurring within 
Ehrenfest time. 

Although the theoretical analysis does not take into 
account either partial transport barriers [40, 41], or the 
’’sticky” dynamics at the regular-chaotic border [ 20 , 21 ], 
which all cause long-time correlations to decay alge¬ 
braically rather than exponentially, we argue in what 
follows that the current model is suitable for the resolu¬ 
tion of our experiment. Figure 5 illustrates the survival 
probability in the chaotic region, obtained from extensive 
ray-dynamics simulations of the microcavity-shaped bil¬ 
liard: despite an overall power-law decay, a closer look at 
the short-time dynamics reveals that the decay is initially 
exponential, behavior known as transient chaos [42]. 

The estimates for the dwelling time and the fitted 
value from the data of 7 , maximum decay rate of the 
chaotic resonances, are both within the timescale of ex¬ 
ponential decay. That suggests that the suppression by 
the absorber of the longest-lived resonances (that scale 
algebraically with 7 [ 21 ]) alone does not affect the num¬ 
ber of excited WGMs measured in the experiment, and 
therefore is not detected by the current apparatus. 
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Supplemental Information to ‘Statistics of chaotic resonances in an optical microcavity’ 

This Supplementary Material is organized as follows. In Sec. I, we describe the boundary shape of the cavity and 
its unidirectional emission. In Sec. II, the Q factor magnitude distribution of the chaotic rays is estimated by ray- 
dynamics simulation. In Sec. Ill, we show geometry and field distribution of microtoroid by SEM images and finite 
element method modeling. In Sec. IV, power loss in the silica waveguide bonding with a thick silicon layer is studied. 
In Sec. V, we provide a cross-check that many high-Q modes do exist in the deformed microcavity by exciting them 
directly through a tapered fiber. In Sec.VI, the process of obtaining the number of high-Q regular modes is described. 
Sec. VII deals with the definition of the Ehrenfest time and rescaling of Lyapunov exponent in the Wey law. In Sec. 
VIII, we explain the test used to assess the goodness of our fits. In Sec. IX, more experimental and fitting results 
are shown. Sec. X contains tables with the quantities and parameters used in our final results. In Sec. XI, the phase 
space with a small deformation factor together with the survival probability of a ray in the chaotic region are shown. 
In Sec. XII, we study the relation between the decay rates of WGMs and their coupling to chaotic modes. 


I. BOUNDARY SHAPE OF THE CAVITY 


The deformed microtoroid cavity in our experiment has a boundary shape given by the curve 


^(A,\ _ / Po(l + e Ei=2.3 COS* (j>) 
^ 1po( 1 + e Ei=2.3 cos* </>) 


for cos 0 > 0 , 
for cos 0 < 0 , 


( 1 ) 


with a 2 = —0.1329, as = 0.0948,^2 = —0.0642, and 63 = —0.0224. The WGMs in the deformed microcavity have 
been demonstrated to possess ultrahigh quality factors in excess of 10 ^ in the 1550 nm wavelength band and to 
exhibit highly directional emission towards the 180° far-held direction, which emits tangentially along the cavity 
boundaries at polar angles 0 = 7 r /2 and 0 = 37 r/ 2 . The deformation is controlled by 77 = (dmax — dmin)/dmax, dmax 
and dmin, respectively, the maximum and minimum diameters of the cavity. The parameter 77 is related to e through 
77 = e |a 2 + as + 62 — 63 ! / 2 . The 77 of cavity shape we used to draw the Poincare surface of section in the main text is 
11.7%. 


H. ESTIMATING THE Q FACTOR MAGNITUDE DISTRIBUTION OF THE CHAOTIC RAYS 
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FIG. 1. Statistics of the quality factor for the chaotic rays in the deformed microcavity, from a ray-dynamics simulation. 

A ray-dynamics simulation is performed on the billiard with boundary set by Eq. 1 (e=l). A swarm of initial 
points, randomly chosen in the chaotic region of the phase space, is iterated 1000 times. Each ray is considered to 
have escaped when it falls below the critical line of total internal reflection, and its quality factor Q is measured as 
Q = 27rrz/, with r lifetime of the ray and v frequency of the light. The statistics shown in Eig. 1 indicates that Q is 
mostly of the order of 10 ^, confirming the estimation from the experimental results. 
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III. GEOMETRY AND FIELD DISTRIBUTION OF THE MICROTOROID 



FIG. 2. (a) SEM cross-section image of the microtoroid. It is an oblique view taken by an angle of 56° with the horizontal 
direction, (b) COMSOL Multiphysics finite element method modeling of a fundamental TE mode. The white solid curve is 
the boundary of the cavity. The scale bar is 3 /xm. 


The field localized at the boundary of the cavity (WGM) is shown in Fig. 2 (toroidal section). It can be seen that 
the toroid has little influence on the propagation of chaotic light fields into the disk part of the cavity since the optical 
WGMs mainly locate at the same plane with the disk. This observation supports our two-dimensional model of the 
micro cavity. 


IV. POWER LOSS IN THE SILICA WAVEGUIDE BONDING WITH A THICK SILICON LAYER 


The finite-element simulation in Fig. 3 shows that the power decrease sharply when the light propagates in the 
2-/im-thick silica waveguide above a silicon layer with a sufficient thickness. In particular, the power decrease by 95% 
when propagating 20 jam (about 19 wavelengths). The wavelength of the light used here is 1550 nm. As for 635 nm 
wavelength, silicon has strong absorption, which may lead to a faster power decrease rate. 




Propagating distance (pm) 


FIG. 3. (a) Finite-element-method modeling the light propagating inside the 2-/xm-thick silica waveguide bonding with a thick 
silicon layer, (b) Remaining energy ratio in silica vs. the propagating distance. 


V. FIBER TAPER AND FREE-SPACE COUPLING 


In Fig. 4(a) we provide an independent cross-check that many high-Q modes {Q > 10^) do exist in the deformed 
microcavity, by exciting them directly through a tapered fiber. Here the silicon pillar attached to the microtoroid 
has largest size, so that dynamical tunneling is inhibited and no WGM can be excited with the free-space coupling 
[Fig. 4(c)]. By reducing the size of the silicon pillar, as shown in Fig. 3 in the main text, high-Q modes can be excited 
indirectly by a free-space beam. 
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FIG. 4. Normalized transmission and top-view images of the cavity coupled by fiber taper [(a) and (b)] and free-space laser 
beam [(c) and (d)]. The absorber-to-cavity ratio r ~ 0.83. 


VI. COUNTING MODES IN TRANSMISSION SPECTRA 

From the transmission spectra collected by the photon receiver, we single out the modes that have relatively high 
Q factors. In particular, we count those with Q > 10^. There are TE modes and TM modes, which have orthogonal 
polarization in the microcavity. A polarization controller is used to excite modes with one polarization, adjust the 
polarization to alternatively achieve the highest efficiencies for one kind of modes and suppress the other one. We 
counted the number of modes in the transmission spectra with both kind of polarizations and add them together 
to obtain the total number of high—Q regular modes. The so called TE modes and TM modes are not perfectly 
orthogonal in the real microcavity, so that some modes may be counted twice, which is the source of the uncertainty 
in our data. 


VII. EHRENFEST TIME AND LYAPUNOV EXPONENT 


We first give a brief explanation of how the semiclassical correction Eq. (4) of the main text was obtained in 
reference [39]. Our first goal is to estimate the fraction of trajectories that escape within the so called Ehrenfest time 
TEhr^ that is the time of quantum-to-classical correspondence. Eor an open system, TEhr is the time it takes for a 
density as large as the size of the channels to be reduced to one Planck cell by the chaotic motion: 


N 

M 


e-f^t ^ 


1 

M’ 


( 2 ) 


verified for TEhr = log V (here fl is proportional to the Lyapunov exponent of the closed system). Equation (2) 
can also be thought of as the mean time it takes for a wavepacket the size of a Planck cell to escape, driven by classical 
dynamics. We notice TEhr is proportional to the logarithm of the number of channels, meaning a faster escape implies 
a longer quantum-to-classical correspondence. Classically, the probability for an initial area Aq of the phase space to 
survive for a time t is 


P(i) ~ (3) 

SO that the probability for an area to survive in the phase space within Ehrenfest time is simply P{TEhr) ^ . 

Now take Aq = M, total area of our phase space in units of h. Erom the definitions we gave of TEhr = log A", 
we can express the survival probability at Ehrenfest time in terms of the number of open channels and Lyapunov 
exponent, as and therefore the number of surviving modes as 

M{TEhr) = (4) 

as seen in Eq. (4) of the main text. In some sense, this is the ‘real’ Weyl law here, since its ‘fractality’ lies in the 
non-integral power. 

The classical estimate of the prefactor involves Ehrenfest time, defined for open systems as 

1 1 th 

TEhr = - log-. 

Td 


( 5 ) 
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Here fi is the Lyapunov exponent of the closed system, Td is the dwelling time that we are already familiar with, while 
th is the Heisenberg time 


th = 


AE' 


(6) 


with AE mean level spacing, that is average distance (or difference) bewteen consecutive energy levels. We know, on 
the other hand, that E = hu, and we may therefore express Heisenberg time in terms of the frequency spacing 


1 


Th = 


Az/’ 


( 7 ) 


and Ehrenfest time as 


^Ehr 


1 N 

ji AT 


(8) 


Here N is the number of open channels as we know, whereas AT = MTAz/, that is the mean frequency spacing times 
the Poincare time (to make it dimensionless), times the number of states M. In plain words, AT is the frequency 
range of our modes in units of the Poincare time. At this point we can still write 


provided that 


TEhr = - log N, 


log N 

log A — log AT ^ 


(9) 


( 10 ) 


Thus we have determined the rescaling to the Lyapunov exponent, following the definition of the Ehrenfest time. 

We may give the rescaling factor an estimate, based on our experimental setup. Taking for example ^ ^ — 
4 • 10“^^s (a is the diameter of the cavity, c the speed of light in the silica), A = 6.35 • 10“^m (visible light), the 
mean wavelength spacing (from the spectra) as AA ^ 10“^^m, we estimate Av = ^ lO^^Hz. Moreover, the total 

number of states is estimated as M Aa/X 40, with A area of the phase space in exam, while the mean number of 
open channels is estimated by the same expression with a smaller area, indicating the open strip of the phase space. 
Using the figures above, we get that 


( 11 ) 

This correction turns out to be within the uncertainties of the average finite-time Lyapunov exponent numerically 
computed for the open billiard, and therefore it has been neglected. 


VIII. TEST 


The expression for the lest used to assess the goodness of our fits is 


X 


2 


1 {Xoh — ^th)^ 


( 12 ) 


Xob being the observed datum, Xth its theoretical expectation, the experimental uncertainty, and e the number of 
degrees of freedom. In general, ^ 1 indicates a poor model fit, while > 1 indicates that the fit has not fully 
captured the data (or that the error variance has been underestimated). In principle, a value of = 1 indicates that 
the extent of the match between observations and estimates is in accord with the error variance. The limit of y^ < 1 
indicates that the model is over-fitting the data: either the model is improperly fitting noise, or the error variance 
has been overestimated. 


IX. MORE EXPERIMENTAL AND FITTING RESULTS 

In Fig. 5 we show more experimental and fitting results. Fig. 5(a)- (c) are in the infrared wavelength band and 
(d),(e) are in the visible wavelength band. All of them support our claim in the main text that the semiclassical 
correction (red solid curves) fits the experimental data better than the purely RMT-based estimate (blue dashed 
curves), especially at smaller deformation, where the two predictions differ the most. 




5 



FIG. 5. (color online) Number of high-Q regular modes (riuj) observed in the transmission spectra of the microcavity (Dots), 
as a function of rescaled absorber-to-cavity ratio Blue dashed and red solid curves are respectively RMT- and semiclassical 
prediction best fits. Here sin^c — 0.69 and sinGh = 0.6. Inset: the area where the two curves differ most. 


X. TABLE OF PARAMETERS 


We report the parameters used in our fittings in the following tables 


TABLE I. Quantities used in relation to the RMT-based expression (Eq.( 6 ) of the main text). 7 is in unit of T with 
T ~ 4 • 10“^^s Poincare time. 


F 

7 

TJ 

A(nm) 

M 


5.2 

0.15 

4 . 2 % 

630 

40 

3.5 

3.2 

0.19 

4 . 2 % 

1550 

20 

3.6 

7.2 

0.16 

6 . 0 % 

630 

40 

5.6 

3.6 

0.18 

6 . 0 % 

1550 

20 

1.8 

4.7 

0.31 

11 . 7 % 

1550 

25 

1.7 


TABLE 11. Quantities used in relation to the semiclassical prediction (Eq.(7) of the main text). Both 7 and are in units of 
T“^, with T ^ 4 • 10“^^s Poincare time. 


F 

7 

r] 

A(nm) 

M 



2.9 

0.19 

4 . 2 % 

630 

40 

0.13 

1.0 

2.0 

0.23 

4 . 2 % 

1550 

20 

0.13 

0.8 

5.0 

0.20 

6 . 0 % 

630 

40 

0.15 

1.0 

2.5 

0.21 

6 . 0 % 

1550 

20 

0.15 

0.8 

1.5 

0.38 

11 . 7 % 

1550 

25 

0.21 

1.3 


XI. PHASE SPACE WITH A SMALL DEFORMATION FACTOR 

The phase space reported in the main text [Fig. 1(b)] belongs to a microcavity with relatively large deformation 
factor, T] = 11.7%. Here we report the classical phase space of the microcavity with the deformation factor r] = A.2% 
(Fig. 6), together with the survival probability of a ray in the chaotic region, obtained by analogously to the result 
shown in Fig. 5 of the main text. One can see here as well an overall algebraic decay of the survival probability, and 
a short-time (transient) hyperbolic behavior. 
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FIG. 6 . (a) Poincare surface of section of the microcavity with deformation factor 77 = 4.2%. Red solid line indicates the angle 
of total internal reflection, (b) (logarithmic scale) Points: average survival probability P(t) of a ray in the microcavity vs. t 
(in units of Poincare time) at ^ = 0.1, from 10® randomly-started trajectories. Line: the ansatz P(t) = exp(—t/r^), Td = 14. 
Inset: the long-time simulation showing algebraic decay. 


XII. DECAY RATES OF WGMS AND COUPLING TO CHAOTIC MODES 


The approximations 


r Yl IX. 

^ In 7 

Y- K? 1/2 

> — ~ riy — . 

V 7n G 


(13) 


leading to Eq. (3) in the main text for the probability of excitation of a WGM may be thought of as too rough. In 
particular, the second expression in (13) means that we are ignoring the fluctuations of both the coupling and the 
decay of the chaotic modes in consideration. In order to show that we are allowed to do so, we shall take a step back 
and rewrite the expression for the amplitude 


Ci^ 


EoEnfn^ 

[lu + - Wo)] + ^ 


Here we recognize the total decay rate of the WGM of frequency uj as 


..tot 


■ luj ■ 


W ~ 'V 

■ > - ~ 7u; ■ 

„ 7n 


1/2 


(14) 


(15) 


where the first term indicates the intrinsic linewidth of the WGM, while the second represents the decay into the 
chaotic modes. That already suggests that the larger n , the larger jtot- If we can really ignore the fluctuations 7 ^ 
in V and 7, we should be able to see that trend. The problem is that we do not count n directly, and thus we need 
to express jtot in terms of some measurable quantity. A good candidate would be then the number of excited 
WGMs. And that reminds us of another important approximation: 


r 2 -h( 


(16) 


where k is unknown. The assumption is that the number of excited WGMs is simply proportional to the probability 
of excitation of one WGM, where, in reality, F = 7 ^ 7 /U^ should be a function of cj, that is even the average coupling 
of each regular mode to the chaotic sea depends on where the mode is supported, in the phase space. It would appear 
from the literature on dynamical tunneling [e.g. A. Backer et ah, Phys. Rev. Lett. 100, 104101 (2008)] that we are not 
allowed to ignore the dependence of F on cj as we did, since the couplings V stretch over several orders of magnitude, 
depending on the regular mode in question. Still, suppose for a moment we can go on making that approximation. 
Eq. (15) would become, in terms of riuj, 


iX = 7c. 







K — n 


UJ 


(17) 
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The advantage of this equation is that we have experimental data to fit it to. Figure 7 shows that Eq. (17) does 
qualitatively capture the behavior of the average linewidths of the WGMs, within some errors. The fitted value 
for 7 ^ corresponds to an average intrinsic Q factor of the order of 10^, which is realistic. Importantly, the overall 
enlargement of the average linewidths with the number of observed WGMs constitutes independent evidence for the 
approximations leading to Eq. (16) to be reasonable for our experiment. 



Number of high-Q modes 

FIG. 7. Points: average linewidth of the excited WGMs vs. their number. Each data point represents one experiment with 
a different size of the silicon pillar. As we know, the number of observed WGMs increases as the size of the pillar decreases. 
Line: Eq. (17), with juj = 326, and k = 171. 





